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Abstract. We have carried out calculations of ionization equilibrium and deuterium fractionation for conditions 
appropriate to a completely depleted, low mass pre-protostellar core, where heavy elements such as C, N, and 
O have vanished from the gas phase and are incorporated in ice mantles frozen on dust grain surfaces. We put 
particular emphasis on the interpretation of recent observations of H2D + towards the centre of the prestellar 
core L 1544 (Caselli et al. 2003) and also compute the ambipolar diffusion timescale. We consider explicitly the 
ortho and para forms of H2, Hjj", and H2D + . Our results show that the ionization degree under such conditions 
depends sensitively on the grain size distribution or, more precisely, on the mean grain surface area per hydrogen 
nucleus. Depending upon this parameter and upon density, the major ion may be H + , Hg~, or D3". We show that 
the abundance of ortho-H2D + observed towards L 1544 can be explained satisfactorily in terms of a complete 
depletion model and that this species is, as a consequence, an important tracer of the kinematics of prestellar 
cores. 
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1. Introduction 

Determining the physical structure of pre-protostellar 
cores is one of the keys to understanding the development 
of protostars. In the simplest view of the development of 
such cores, they contract slowly towards some "pivotal 
state", subsequent to which dynamical collapse sets in. 
Clearly, one should attempt to determine the parameters 
of this pivotal state, since these parameters dictate the 
subsequent evolution. It seems probable that the pivotal 
state is marked by high column density and low temper- 
ature. Likely candidates can be identified using measure- 
ments of dust emission and absorption (e.g. Andre et al. 
2000), which indicate temperatures (of both gas and dust) 
of about 10 K or below, molecular hydrogen column den- 
sities in the range 10 22 to 10 23 cm" 2 , and sizes of several 
thousand AU. 

An additional result from recent observational studies 
of prestellar cores is that depletion of molecular species 
on to dust grain surfaces is also a marker of relatively 
evolved cores (i.e. close to the pivotal state) with high 
central densities and column densities (e.g. Caselli et al. 
2002a, Tafalla et al. 2002). Tafalla et al. show that CO 
becomes depleted by at least a factor of 10, relative to 
its canonical abundance, above densities of 5 10 4 cm" 3 . 
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CS shows similar behaviour, whereas nitrogen containing 
species such as NH3 and N2H+ have abundances which 
are either constant or (in the case of ammonia) increase 
somewhat in the high density gas surrounding the dust 
emission peak. Their interpretation of these observations 
is that molecular nitrogen (N2), which is the source of 
the observed NH3 and N2H + , is sufficiently volatile to re- 
main in the gas phase at densities of order 10 5 cm -3 (es- 
sentially because of spot heating of the grain surfaces by 
cosmic rays) , whereas CO (the source of carbon for most 
C-containing species) is not. 

What happens at still higher densities? It seems likely 
that, at sufficiently high densities, all species containing 
heavy elements will condense out. The difference in sub- 
limation energies for CO and N2 is not great (e.g. Bergin 
and Langer 1997) and thus one might expect N2 and other 
N-containing species to disappear from the gas phase 
at densities only somewhat higher (or dust temperature 
slightly lower) than found for CO. Whilst it is true that 
the laboratory experiments may not accurately simulate 
interstellar ice surfaces, it seems probable that N 2 disap- 
pears from the gas phase at densities above 10 6 cm" 3 . 

In fact, the densities inferred from mm dust continuum 
emission in several objects are of order 10 6 cm" 3 and there 
is some evidence for "holes" in the N2H+ distribution at 
densities above 3 10 5 cm" 3 . For example, Bergin et al. 
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(2002) find a flattening of the N 2 H+ column density dis- 
tribution in B68 which could be interpreted as being due 
to an absence of N 2 towards the dust peak. Bcllochc (2002) 
finds a "N 2 H + hole" around the peak of the Class source 
IRAM04191, which is most easily interpreted in terms of 
N 2 condensing out at the highest densities in the prestellar 
core. Last, but not least, Caselli et al. (2003, in what fol- 
lows CvTCB) have detected compact emission (2000 AU 
in size) in the lio — In H 2 D + line towards the dust emis- 
sion peak of L 1544. This latter discovery is significant in 
view of the fact that CvTCB infer an abundance of order 
10~ 9 relative to H 2 . Estimates of the ionization degree in 
L 1544 (Caselli et al. 2002a) are of order 2 10~ 9 . Taken to- 
gether, these results imply that H 2 D+(and hence probably 
Hg") is a major ion. It seems likely that this can be the 
case only if ions containing heavy elements, such as N 2 H + 
and HCO + are absent, which in turn can be true only if 
the heavy element content of the gas has condensed on to 
grain surfaces. If this is the case, one must rely on species 
lacking heavy elements (such as H 2 D+) to trace the kine- 
matics of the high density nucleus of the pre-protostellar 
core. Also, it becomes relevant to establish the ionization 
degree and other physical parameters of a region where 
species containing heavy elements have condensed out. 

In this paper, we study some of the consequences of 
complete depletion in an isolated, low mass prestellar core. 
In particular, we compute the ionization degree and its de- 
pendence upon physical parameters such as the gas den- 
sity and the grain size. We consider in detail deuterium 
fractionation in these conditions and try to explain the ob- 
served H 2 D + (or, more precisely, the ortho-H 2 D + ) abun- 
dance in objects such as L 1544. In section 2, we discuss 
the model assumptions and, in section 3, we summarize 
the results. Then, in section 4, we consider some of the im- 
plications of our calculations and in section 5, we briefly 
summarize our conclusions. Some preliminary results of 
our study have been presented elsewhere (Walmsley et al. 
2004). 

2. Model 

Our model derives from the work of Flower et al. (2003) 
and Flower and Pineau dcs Forets (2003). The time- 
dependent chemical rate equations are solved numerically 
for a situation in which heavy elements have condensed 
out on to dust grain surfaces. As in the case of the primor- 
dial gas, the only significant gas-phase species are com- 
pounds containing the elements H, D and He. An impor- 
tant role is played by the grains which influence the charge 
distribution (ions, charged grains, and electrons) and en- 
able the formation of H 2 and HD. Apart from these latter 
two processes, we neglect surface chemistry in the present 
study. We discuss below our assumptions concerning re- 
action rate coefficients. 

First, we consider briefly what is implied by the state- 
ment that "no heavy elements are present in the gas 
phase" . For the purpose of the present study, this means 
in essence that ions such as H^, produced as a conse- 



quence of cosmic ray ionization of molecular hydrogen, 
are destroyed by dissociative recombination or recombi- 
nation on grain surfaces, rather than by transferring a 
proton to species such as CO and N 2 which have larger 
proton affinities than H 2 . Assuming a rate coefficient k e for 
dissociative recombination of Hj with electrons of 4 10 -7 
cm 3 s" 1 (McCall et al. 2003, extrapolated to 10 K with 
a T~ - 52 dependence: see Appendix A) and a rate coeffi- 
cient for proton transfer reactions between Hjj~ and species 
X such as CO or N 2 of the order of the Langevin value, 
«i = 10~ 9 cm 3 s _1 , the condition for "complete deple- 
tion" is [X] < fc e [e]/a!L ; where the square brackets de- 
note a fractional abundance; this implies that, for exam- 
ple, the CO abundance [CO] in the "complete depletion 
limit" should be less than 400[c], or roughly 4 10~ 6 in cores 
such as L 1544, where [e] w 10~ 8 . Thus, the calculations 
reported here will be relevant when species such as CO 
have abundances much less than 10~ 6 . 

In what follows, we report the steady-state conditions 
predicted by our calculations. The timescales to reach 
chemical equilibrium are somewhat shorter than the dy- 
namical timescales in cores such as L 1544. The free-fall 
time is approximately 3 10 4 y for a density of 10 6 cm -3 ; 
but, when magnetic fields are present, the collapse time is 
about an order of magnitude larger (see Ciolek and Basu 
2000, Aikawa et al. 2003). As shown below, the equilib- 
rium timescale for the entire chemistry is roughly 10 5 y, 
although ionization equilibrium (with a characteristic time 
l/[k e n(e)], which is of the order of a few years) is reached 
much more rapidly. We conclude that the steady-state as- 
sumption is a valid first approximation but that some of 
our results may need to be reconsidered in the light of 
their time dependence. 

2.1. Grain processes 

In general, we have followed the approach and adopted 
the rate coefficients relating to "large grains" of Flower 
and Pineau des Forets (2003, in what follows FPdesF; 
see their Appendix A). The following processes affecting 
the grain charge are taken into account: electron attach- 
ment; electron detachment by the H 2 fluorescence pho- 
tons, generated by collisions with the secondary electrons 
produced by cosmic ray ionization of hydrogen; attach- 
ment of positive ions to neutral grains and recombination 
of positive ions with negatively charged grains. However, 
for the present applications, their assumption of a Mathis 
et al. (1977, hereafter MRN) power law distribution of 
grain sizes (based on the extinction curve derived for dif- 
fuse interstellar gas) seemed inappropriate. We have there- 
fore assumed single-size grains, with a dust-to-gas mass 
density ratio (0.013) consistent with all elements heavier 
than helium having condensed out. Moreover, we adopt 
a grain material density of 2 gm cm~ 3 for the core plus 
"dirty ice" mantle. We treat the grain radius a g as a pa- 
rameter of the model. 
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Recombination of ions with electrons on negatively 
charged grain surfaces is important because, if the 
available grain surface area is large enough, this pro- 
cess can dominate the destruction of ions such as H^". 
Recombination on grains is certainly the main destruction 
mechanism for protons, which recombine only slowly (ra- 
diatively) with electrons in the gas phase, and hence our 
results depend on the surface area of negatively charged 
grains. Accordingly, we have considered carefully the dis- 
tribution of grain charge, adopting the rate coefficients 
of FPdesF. We assumed that H2 and HD form on grain 
surfaces, at rates which scale with the total grain surface 
area. 

2.2. Chemical scheme and D-fractionation 

The assumption of complete depletion of the heavy ele- 
ments on to grains leads to a considerable reduction in 
the number of chemical species which have to be taken 
into account: H, H 2 , H + , H 2 , H;j~ (and their deuterated 
forms), He and He + . At the low temperatures which pre- 
vail in protostellar cores, deuterium fractionation occurs, 
resulting in greatly enhanced abundances of H2D + , D2H+ 
and D^" . In this situation, the rates of the associated for- 
wards and reverse deuterium-substitution reactions de- 
pend not only on the differences in zero-point energies 
along the deuteration sequence (Ramanlal et al. 2003) but 
also on the ortho/para abundance ratios of H2, H^ , H^", 
and H 2 D + . Accordingly, we have treated the ortho and 
para forms of these molecules as separate species in our 
calculations. 

A key role is played by the ortho/para ratio of molec- 
ular hydrogen, which is expected to be much higher than 
in LTE (see the discussion of Le Bourlot 1991); this has 
consequences for the ortho/para ratios of species such as 
H 2 D+ (see Gerlich et al. 2002, hereafter GHR; Pagani et 
al. 1992), which become enhanced roughly in proportion 
to that of H2. Also important in this context are the spin 
selection rules for ortho-para conversion, discussed, for 
example, by Uy et al. (1997). 

Caution has to be exercised when using the data of 
Uy et al. (1997). These authors considered conditions in 
which the ortho and para forms are statistically popu- 
lated, i.e. in proportion to their nuclear spin statistical 
weights, (21 +1). This situation obtains at temperatures 
and densities much higher than those considered here. In 
prestellar cores, only the low-lying states of the ortho and 
para species are significantly populated. Under these cir- 
cumstances, the statistical weights of the rotational levels, 
(2 J + 1), must be taken explicitly into account. When 
compiling the rate coefficients for reactions involving the 
ortho and para forms of H 2 , H+, H+, and H 2 D+, we made 
the simplifying assumption that only the lowest rotational 
level is populated in each case and treated the correspond- 
ing ortho and para forms as distinct species, which inter- 
convert through chemical reactions. The complete set of 
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reactions included in the model and their rate coefficients 
are given in Appendix A. 

In the limit of complete depiction, deuterium fraction- 
ation occurs mainly in the reactions 

H+ +HD -> H 2 D+ +H 2 (1) 
with rate coefficient ki, 

H 2 D+ +HD -> D 2 H+ + H 2 (2) 
with rate coefficient fc 2 , and 

D 2 H+ + HD — > + H 2 (3) 

with rate coefficient k%. These reactions are exothermic by 
roughly 200 K and consequently enhance abundance ratios 
such as n(H 2 D + )/7i(H;}~). However, such enhancement can 
be limited by a variety of processes. If heavy elements are 
present in the gas phase, the high abundances of molecules 
with high proton affinities, such as CO, result in the frac- 
tionation being transferred to ions such as DCO + ; this 
limits the degree of enhancement of deuterium in H 2 D + . 
It follows that the depletion of heavy elements favours deu- 
terium fractionation (see Dalgarno and Lepp 1984), and 
so a "completely depleted core" is a limiting case where 
one might expect high degrees of enhancement. 

Even in the absence of heavy elements, there are 
several processes which limit deuterium substitution in 
species such as H3. For example, it has been known for a 
long time that dissociative recombination with electrons, 
at a rate fc e n(e), competes with deuterium enrichment in 
H2D" 1 ", for example (Watson 1978). It is also clear that re- 
combination of H2D+, D2H+ and D^~ on grain surfaces can 
have a similar effect, and so models with a large grain sur- 
face area tend to have moderate D-fractionation. The rela- 
tive importance of these two processes is roughly the same 
for D-fractionation and for ion neutralization: in models 
where the D-fractionation is limited by recombination on 
grain surfaces, the degree of ionization is limited by the 
same process. 

D-fractionation in an ion such as H2D + is determined 
also by the transfer of the deuterium to the higher mem- 
bers of the chain, D2H + and Djj\ Thus reactions [3 and 01 
restrict the fractionation in H2D+ and D2H + , respectively. 
As we shall see below, the effect of these reactions is to 
limit ratios such as [H 2 D+]/[H+] and [D 2 H+]/[H2D+] to 
values of the order of unity. 

GHR pointed out that ortho-H2 can restrict deu- 
terium fractionation because the reaction between ortho- 
H2 and ortho-H2D + is exothermic and believed to be fast. 
Depending on internal excitation, this process may also be 
important for D 2 H + . We find that the reaction with ortho- 
H2 is not usually dominant but is significant because the 
abundance of ortho-H2 is high. 

The ratio of the H2D + and H^" abundances is given 
approximately by 

[H 2 D+]/[H+] = (4) 

12 J/L 3J fe e [e] + fc ig [g-]+2fc 2 [D]+2fcr W 
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Fig. 1. Sample result for our reference model (n(H 2 ) = 
10 6 cm" 3 , T = 10 K, a g = 0.1 /xm, C = 310~ 17 s" 1 ) 
showing the approach to steady-state conditions for or- 
tho/para H2, D and H. 

where [D] is the fractional abundance of D in the gas 
phase, ki g is the rate coefficient for neutralization of 
H 2 D + on negatively charged grains of number density 
n g- — [g ] n n, and k^ is the rate coefficient for the re- 
verse of reaction ^ We see that deuterium fractionation 
is limited by the terms in the denominator of the right 
hand side of this equation. One may note also that, whilst 
[H 2 D+]/[H+] and [D 2 H+]/[H 2 D+] are limited to k 1 /k 2 
and kijkj 11 respectively, there is no such restriction for 
p3~]/[D2H + ] because D^" does not react with HD. Thus, 
in certain circumstances, Djj" can be the principal ion. 



2.3. Numerical method 

The code described by Flower et al. (2003) has been used 
to calculate the steady-state abundances of the chemical 
species, Xj. The program integrates the coupled differen- 
tial equations 



dt 



n(Xi) = /[npg.T.t] 



for the number densities n(Xi), as functions of time, t, 
until steady state is attained. The kinetic temperature, 
T, was assumed constant. By adopting this approach, not 
only are the steady-state values of the number densities 
obtained for sufficiently large t but also information is re- 
covered on the timescales required for the various chemical 
reactions to reach equilibrium. 

In Fig. ^ we show the results of a sample calcula- 
tion in which, initially, the hydrogen and deuterium are 
in molecular form, with an elemental abundance ratio 
n-o/nn — 1.6 x 10~ 5 . The initial ortho/para H 2 ratio was 
3.5 x 10~ 7 , the equilibrium value at the gas kinetic tem- 
perature T = 10 K (assumed constant). We find that the 
time required to reach steady state for all species is of the 
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Fig. 2. Abundances of major ions and electrons in our 
standard model (a g = 0.1 fim, T = 10 K, £ = 3 x 
10~ 17 s~ 1 ) as functions of the density of molecular hydro- 
gen. 

order of 10 5 y (almost independent of number density). 
However, as noted above, the time to attain ionization 
equilibrium (a few years) is much shorter than that re- 
quired to attain equilibrium between ortho- and para-H 2 
(the longest chemical timescale). Atomic hydrogen and 
deuterium are intermediate, reflecting the timescale for 
the formation of H 2 and HD on grains. 

The time characterizing chemical steady state is com- 
parable with the free-fall time, rg, at the lowest density 
(riH = 2 x 10 5 cm -3 ) and is an order of magnitude larger 
than Tff at the highest density (riH — 2 x 10 7 cm -3 ) con- 
sidered. However, magnetic support to the cloud increases 
the time for core formation and, in practice, the timescale 
for chemical equilibrium is likely to be comparable with 
that for contraction and collapse of the protostellar core. 

3. Results 

3.1. Dependence on model parameters 

We present our results with reference to a model with the 
following parameters : a g — 0.1 /im, p g = 2 g cm~ 3 , 
T = 10 K, and C = 3 x lO" 17 s -1 , where a g is the 
grain radius, p g is the mean density of the grain mate- 
rial (core and ice mantle) , T is the gas temperature, and C 
is the rate of cosmic ray ionisation of H 2 . To these parame- 
ters correspond a fractional abundance of grains n g /nn = 
3.5 x 10~ 12 and a grain surface area per hydrogen nucleus 
n g a g /nn = 1.1 10~ 21 cm 2 . These numbers have been de- 
rived assuming that all heavy elements are incorporated 
into either the cores or the ice mantles of the grains (cor- 
responding to a dust-to-gas mass density ratio of 0.013). 

In Fig. [21 are plotted the fractional abundances of 
the principal ions and of the electrons as functions of 
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Fig. 3. Ortho/para ratios of H 2 (multiplied by 10 4 ), Hj", , , ... n 

i tt r • r ■ r r \ f model with a smaller gram size ol 0.025 /im. 

and H 2 D + as functions of density lor our reference model 

(cf. Fig. |2J . The approximate expression of Gerlich et al. 

(2002) for H 2 D + is shown for comparison. 



T1QI2) ~ ^h/2; the results of computations of ortho/para 
ratios for the same model are presented in Fig. [3] and will 
be discussed in section 3.4 below. At the lowest density 
considered, n(H 2 ) = fO 5 cm~ 3 , H + is the major ion and 
the removal of molecular ions by dissociative recombina- 
tion is significant. H + , on the other hand, is removed by 
recombination on grain surfaces. As n(H 2 ) increases, re- 
combination on grains dominates dissociative recombina- 
tion for the molecular ions. H + becomes a relatively minor 
species at high density. 

The balance between the various deuterated forms of 
is determined approximately by equation 0] and its 
analogues for [D 2 H+]/[H 2 D+] and [D+]/[D 2 H+]. At densi- 
ties below 10 6 cm -3 , the dissociative recombination term 
k e [e] is dominant in the denominator of the right hand 
side of equation 01 whereas, at higher densities, the term 
due to the reaction with HD becomes more important and 
results in [H 2 D+]/[H+] (and [D 2 H+]/[H 2 D+]) being close 
to unity. On the other hand, D^" cannot be destroyed by 
HD and becomes the dominant ion. 

Fig. is analogous to Fig.0 for an assumed grain ra- 
dius of 0.025 /im. In this case, recombination of ions on 
grain surfaces is the dominant neutralization process for 
all ions and H + is never dominant. Moreover, for densities 
below 10 7 cm -3 , H^ is more abundant than its deuterated 
forms. We note that the increased rate of grain neutral- 
ization causes the ion abundance to fall and its density 
dependence to become steeper. 




3.2. Ionization degree for complete depletion 

The determination of the ionization degree in the region of 
complete depletion is important for a number of reasons. 



n(H 2 ) (cm J ) 

Fig. 5. The free electron abundance computed as a func- 
tion of density for grain sizes between 0.025 and 0.2 ^m. 
The result of McKee (1989) is shown for comparison (bro- 
ken curve). Note the decrease in the level of ionization for 
small grain sizes (large grain surface areas). 

The degree of ionization determines the coupling of the 
magnetic field to the gas and hence the timescale for am- 
bipolar diffusion. The electron abundance plays a crucial 
role also in the chemistry and in deuterium fractionation. 

In Fig. |SJ the free electron abundance, n(e)/n(H 2 ), 
computed as a function of density, is compared with 
the relationship of McKee (1989), n(e)/n(H 2 ) = 1.3 x 
10 -5 n(i/ 2 )~ - 5 , which is often used in discussions of core 
evolution. It may be seen that our model generally pre- 
dicts lower degrees of ionization. We find that n(e)/n(H 2 ) 
varies approximately as n^ ' 75 between rin = 2 x 10 6 and 
2 x 10 7 cm~ 3 . As already noted, the degree of ionization 
of the gas is sensitive to the assumed grain size, owing 
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10 5 10 6 10 7 

n(H 2 ) (cm 3 ) 

Fig. 6. Ratio of ambipolar diffusion time tad (equation^ 
to free fall time as function of molecular hydrogen den- 
sity for values of the grain radius between 0.025 and 0.2 

[ITCl. 

to the importance of the recombination of ions with elec- 
trons on grain surfaces. The fraction of electrons attached 
to grains remains small for all the grain sizes considered 
here: the ratio of negatively charged grains to free elec- 
trons attains 0.20 for for a g — 0.025 /xm and n(H2) = 10 7 
cm -3 . However, we note that this ratio becomes larger for 
smaller grain sizes, becoming >> 1 for a g < 0.01 /im. The 
fraction of grains which is neutral decreases from 0.68 for 
a g = 0.025 ^m to 0.20 for a g = 0.2 ^m and for n(H 2 ) = 10 7 
cm~ 3 (but this ratio is insensitive to the gas density). The 
fraction of positively charged grains is always negligible. 

3.3. Ambipolar Diffusion in depleted cores 

If the core is sub-critical (i.e. if the magnetic field is capa- 
ble of resisting collapse), the timescale for star formation 
is determined by that for ambipolar diffusion. Here we 
estimate the timescale for ambipolar diffusion under the 
conditions of complete heavy element depletion. 

The ambipolar diffusion timescale, tad, depends on 
the ion abundance, rii/n n , the rate coefficient for momen- 
tum transfer between the ions and the neutrals (essentially 
H2), < civ >m~ 2 x 10 -9 cm 3 s -1 , and the masses m» and 
m n of the charged and neutral collision partners. We have 
taken 

2 \ rii m,im n 

tad « n 2 > : < av > m (5) 

7rG m„ *■ — ' n n m* + m n 

n i 

where G is the gravitational constant and the summation 
is over the ionic species. The rate coefficient for momen- 
tum transfer approaches its Langevin form 

< (TV >i n = 2-Ke(a/m in ) - 5 



at low temperatures (cf. Flower 2000), and hence tad 
becomes proportional to the square root of the reduced 
mass, rrii n — [m,im n / (m, + m n )]; a — 7.7 x 10 -25 cm 3 is 
the polarizability of H2. Thus, tad increases by 50 per- 
cent as the dominant ion changes from being H + to being 
D^" . In Fig. El we plot the ratio of the ambipolar diffusion 
timescale tad to the free-fall timescale t// as a function 
of density. We see from Fig. that the ambipolar diffusion 
timescale in the limit of complete depletion is of the same 
order of magnitude as that inferred for normal ionic abun- 
dances in dense cores, i.e. about an order of magnitude 
larger than the free-fall time (see, for example, Shu et al. 
1987, Ciolek and Basu 2000). As one might expect, smaller 
grain sizes lead to shorter ambipolar diffusion timescales, 
owing to more rapid neutralization on grain surfaces. 

The above discussion neglects the collisional coupling 
between the neutral gas and charged grains (cf. Nakano 
and Umebayashi 1980). For the conditions that we have 
considered, the coupling to the charged grains is domi- 
nant, particularly for small grain sizes and high gas den- 
sities. We estimate that its effect is to increase the value 
of tad by one to two orders of magnitude. We conclude 
that the collapse of a magnetically sub-critical, completely 
depleted core is likely to be impossible on reasonable 
timescales. 

3.4. The abundances of ortho-ty, ortho-fyD+and 
ortho-H^ 

Fig. |3 shows that the ortho/para ratios considered here 
are relatively insensitive to density. For molecular hydro- 
gen, we derive a steady-state ortho/para ratio of approxi- 
mately 5 10 ~ 5 for the conditions in the nucleus of L 1544, 
and this value is reflected in the H2D" 1 " ortho/para ra- 
tio, as shown by GHR. In fact, we find that the H2D 4 " 
ortho/para ratio predicted by equation (7) of GHR is a 
good approximation to our numerical results. 

The expression of Le Bourlot (1991) for ortho/para 
H2 was derived for different conditions (lower density and 
high abundances of molecular ions such as HCO + ) and 
is not a good approximation in the conditions that we 
have considered. In our situation, the ortho/para H2 ratio 
is given approximately by C/[^op w(XH )], where n(XH + ) 
represents the sum of the number densities of ions such 
as H + and H^ which convert ortho- to para-H2 through 
proton transfer reactions, with a rate coefficient for ortho 
to para conversion k op which is of order 10 -10 cm -3 s — 1 . 

The ortho/para H2D 4 " ratio has a flat dependence on 
density, reflecting the behaviour of ortho/para H2. The 
ortho forms of both H2D + and H3" are produced mainly 
in reactions of the para forms with ortho-H2 and so their 
high ortho/para ratios are attributable to the relatively 
high ortho-H2 abundance; Fig. illustrates this point. 
Because the ortho/para H2 ratio is not thermalized under 
the conditions which we consider, the ortho/para H2D 4 " 
and H^" ratios are not thermalized either. It follows that 
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Fig. 7. Ortho/para ratios of and H2D" 1 " as functions of 
the molecular hydrogen ortho/para ratio. The conditions 
are those of our reference model: n.(H2)= 10 6 cm~ 3 , T = 
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Fig. 8. The computed ortho-H2D + abundance from our 
reference model (T = 10 K and ( = 3 x 10~ 17 s _1 ) as a 
function of density for several values of the assumed grain 
radius. The values observed by Caselli et al. (2003) are 
shown for comparison. 



measurements of the ortho/para H3 ratio should not be 
used to infer a temperature. 

The basic question that we wish to answer is whether 
the ortho-H2D + emission observed in L 1544 can be re- 
produced by a model which assumes complete heavy el- 
ement depletion. The abundance of ortho-H2D + relative 
to H2 observed towards the peak of L 1544 is 5 10 -10 
according to CvTCB. In Fig. |H| are shown the computed 
values of [ortho-H2D + ]/[H2] as a function of density for 
various assumed grain sizes. There is a tendency for [ortho- 
tbD+J/prb] to decrease with increasing density, owing to 
the fact that deuteration proceeds further down the chain 
(towards D3 ) at high density. Our results are within a 
factor of 2 of the observed peak ortho-H^D" 1 " abundance, 
for grain radii below 0.1 /mi. Given the observational un- 
certainties, this level of agreement might be considered 
satisfactory. It is also significant that computations for 
a = 0.01 /im show that the computed ortho-H2D + abun- 
dance decreases for all values of n(H2) relative to the re- 
sults in Fig. |H1 for a = 0.025 /im, and thus grain radii in 
the range 0.025-0.05 /im maximize the ortho-H2D + abun- 
dance. We note that the "observed" abundance assumes 
a rate of collisional excitation of H2D + , a temperature, 
and a density distribution, of which the latter two are 
controversial (see, for example, Galli et al. 2002). Higher 
angular resolution studies of the mm-submm continuum 
are needed for further progress. 

One result from our study is that a relatively high rate 
of cosmic ray ionization is required in the nucleus of L 1544 
to account for the observed ortho-H2D + emission; this im- 
plies the penetration of low energy (100 MeV) cosmic rays 
into the high density core. Cosmic ray fluxes higher than 
assumed here would enhance the computed ortho-H2D + 
abundance; but to a higher flux would correspond an in- 



creased rate of heat input to the gas and temperatures in 
excess of those observed in NH3 and other tracers. 

3.5. Atomic Deuterium in the completely depleted core 

An interesting by-product of our work is the prediction 
of the atomic D abundance under the conditions of com- 
plete depletion. Fig. [5] shows the abundances of atomic 
deuterium and hydrogen as functions of density for vari- 
ous grain radii. Atomic H is produced mainly by cosmic 
ray ionization and destroyed by the formation of H2 on 
grains. As a result, the fractional abundance of H varies 
inversely with the density and decreases with the grain ra- 
dius, owing to the greater grain surface area available to 
form H2. Atomic D, on the other hand, is produced by the 
recombination of dcuterated ions and becomes relatively 
more abundant as deuterium fractionation increases. It is 
noteworthy that the abundance of atomic D found from 
our computations [> 0.1 cm~ 3 for n(H 2 ) = 10 6 cm -3 ] is 
of the order of that required to explain the high fractiona- 
tion of deuterium in methanol and formaldehyde on grain 
surfaces (Caselli et al. 2002c). 

4. Discussion 

4.1. Grain size distribution in depleted cores 

The grain size distribution in prestellar cores is not known. 
It depends on the timescales for grain coagulation and 
mantle accretion, relative to the dynamical and chemi- 
cal timescales. In view of this uncertainty, we compared 
the results of computations assuming a unique grain ra- 
dius with the corresponding calculations assuming a MRN 
power-law size distribution; we retained the same value 
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Fig. 9. Atomic H (full lines) and atomic D (broken lines) 
densities for four values of the grain radius (increasing 
from bottom to top) for densities between 10 5 and 10 7 
cm~ 3 . 

of the grain surface area per H-nucleus, n g a g /nn, in 
both cases. In calculations with a unique value of a g , 
n g <Tg/nn = 1.1 1CT 21 (0.10/a s (/im)) cm 2 . When a power 
law dependence a~ r is adopted, providing the exponent 
r > 3 (MRN derived r = 3.5), the lower limit to the ra- 
dius essentially determines the grain surface area. When 
r = 3.5, the lower limit to a g is 600 A for n g a g /na = 
1.1 10~ 21 cm 2 . We found that the results of calculations 
with the same value of grain surface area per H-nucleus 
were quantitatively similar for the cases of a unique grain 
radius and a MRN power-law size distribution. 

As shown above, the detection of H2D + in L 1544 
by CvTCB can be explained if a g < 0.1 /im or, equiv- 
alently, if the grain surface area per H-nucleus is at least 
1.1 10~ 21 cm 2 , which is close to the value for grains in 
diffuse clouds. This conclusion is consistent with deduc- 
tions which have been made from the mm-infrared extinc- 
tion curve in prestellar regions: Bianchi et al. (2003) found 
that the ratio of mm to NIR dust absorption is marginally 
larger in prestellar environments than in diffuse clouds; it 
is also consistent with the simulations of Ossenkopf and 
Henning (1994). 

Under conditions of complete heavy element depletion, 
the thickness of the grain mantle is expected to be roughly 
40 percent of the grain radius, or 400 A for 0.1 /im grains. 
Then, 45 percent of the grain mass is in the refractory core 
and 55 percent in the ice mantle. A MRN size distribution 
yields an ice mantle 260 A thick, independent of the grain 
radius: see Appendix B. 

4.2. Ortho to para ratios 

Our calculations show that the ortho/para ratios of both 
Hjj" and H2D+ are directly dependent on the H2 or- 



tho/para ratio. Thus, measurements of the former would 
yield information on the latter. We emphasize that, whilst 
steady state may be a reasonable approximation for most 
of the chemical species that we have discussed, this may 
not be the case for the H2 ortho/para ratio. Thus, reli- 
able measurements of the ortho/para ratios of either H^" 
or H2D + could provide information on the evolution of the 
source. 

Unfortunately, the loi — 0oo transition to the ground 
state of para-H2D + (at 1.37 THz) will be difficult to ex- 
cite at the temperatures of prestellar cores and may be 
observable only in absorption towards a background con- 
tinuum source. On the other hand, Hjj" has no allowed 
submm transitions and is best observed in the NIR at 3.5 
/im in absorption towards background sources (McCall et 
al. 1999). It has been assumed that the H^" which has been 
observed exists in regions whose composition is similar to 
that of "normal" interstellar molecular clouds and, fur- 
thermore, that the ortho/para forms are in LTE. In our 
view, both of these assumptions may prove to be invalid. 
Observations of H2D+ and the other members of the fam- 
ily may help to shed light on this issue. 

4.3. Observations of completely depleted cores 

Observations of pre-protostellar cores, in which heavy el- 
ements have frozen on to grains, are important because 
they provide us with information on the structure and 
kinematics of the core at the time of its initial gravita- 
tional collapse. It is important to observe the mm contin- 
uum emission of these objects with high angular resolution 
in order to derive the density and temperature distribu- 
tions for comparison with theoretical determinations (e.g. 
Galli et al. 2003). Likewise, kinematical information may 
be obtained from higher angular resolution studies of the 
lio — In hne of H 2 D+. The present study shows that 
it is plausible to conclude that the ortho-H2D + emission 
detected by CvTCB towards L 1544 arises in completely 
depleted material. Indeed, these authors failed by a factor 
of 10 to account for their observations of L 1544 using es- 
timates based on less extreme assumptions regarding the 
depletion of heavy elements from the gas phase. 

Our results can be compared with those of Roberts et 
al. (2003) who carried out a time dependent calculation of 
the chemistry in a prestellar core with density 3 10 6 cm~ 3 
including freeze-out. They found, as do we, that can 
be the most abundant ion, though their results suggest 
that H^" becomes the major ion when all heavy species 
freeze out; this may be attributable to a different way of 
dealing with grain neutralization of ions. 

Further progress in this area will depend on our ca- 
pacity to sample transitions other than the the lio — In 
line of H2D" 1 ". H^ and D;j" do not have permanent dipole 
moments, and hence their detection is feasible only in ab- 
sorption in the near infrared. Our calculations suggest that 
searches for D^" should be directed towards regions where 
high heavy element depletion is believed to have occurred. 
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Searches for emission in low-lying transitions of D 2 H + 
should also be made: the abundance of D 2 H+ in condi- 
tions of complete depletion is predicted to be comparable 
to that of H 2 D+. Under conditions of LTE at temperature 
T, the relative populations of the lowest ortho (0 o) and 
para (loi) levels of D 2 H + would be 

= (9/6)exp(-50.2/T) 

where the nuclear spin statistical weights (21 + 1) of the 
ortho (I = 0, 2) and the para (I = 1) levels are 6 and 
3 and their rotational degeneracies (2 J + 1) are 1 and 3, 
respectively, yielding total statistical weights of 6 and 9 for 
Ooo and loi, respectively. However, under the conditions 
in prestellar cores, the para/ortho D 2 H + ratio (like the 
ortho/para ratios discussed above, and for related reasons) 
is probably much in excess of its LTE value; we estimate 
that the para and ortho forms of this species may have 
comparable abundances. Thus a deep search for the 692 
GHz lio — loi transition of D 2 H + seems worthwhile. 

5. Conclusions 

We summarize here our most important conclusions. 

— We believe that our data strongly suggest that the de- 
tection of H 2 D + by CvTCB is due to complete freeze- 
out of the heavy elements in the core of L 1544. If so, 
only species such as Hjj" and its deuterated counter- 
parts are useful tracers of the kinematics in the dense 
centre of L 1544 and other similar objects. 

— Knowledge of the mean grain surface area is needed 
in order to understand phenomena such as deuterium 
fractionation; it seems likely that this remains true 
even if some of the heavy species are in the gas phase. 
Our results imply that, in L 1544, the mean grain sur- 
face area per H-atom is above 1CP 21 cm 2 , which is 
comparable to values found in the diffuse interstellar 
medium. 

— The degree of ionization and the timescale for ambipo- 
lar diffusion, determined under conditions of complete 
heavy element depletion, are not significantly different 
from values reported elsewhere in the literature, eval- 
uated assuming degrees of depletion more typical of 
dense interstellar clouds. However, we estimate that, 
under the conditions considered here, the effect of the 
coupling between the neutral gas and charged grains 
is to increase the ambipolar diffusion timescale by be- 
tween one and two orders of magnitude; this finding 
will be the subject of further investigation. 

Appendix A: Chemical reactions and rate 
coefficients 

In this Appendix are specified the reactions included and 
the rate coefficients adopted in our chemical model. 

Protons have nuclear spin 1=1/2 and are fermions. 
The ortho and para forms of H 2 , Hj , Hg and H 2 D+ have 



been treated as distinct species. Ortho- and para-H 2 were 
assumed to form hot, on grains, in the ratio 3:1 of their nu- 
clear spin statistical weights. Subsequent proton-exchange 
reactions interconvert the ortho and para forms. Such re- 
actions were assumed to involve only the rotational ground 
states of the ortho and para modifications, and the cor- 
responding rate coefficients take account of the relative 
statistical weights (27+1) (2 J+l) of the products and re- 
act ants. 

Deuteration of H J in reactions with HD may be viewed 
as proton transfer (from Hg to HD) or dcuteron-proton 
exchange; we allowed for both viewpoints in our reaction 
set. Rate coefficients for dissociative recombination, with 
electrons in the gas phase or on negatively charged grains, 
were assumed to be the same for the ortho and para forms 
of the recombining ion. 

We note that deuterons have nuclear spin 1=1 and 
are bosons. Multiply deuterated species have ortho and 
para forms. In the present study, we have not attempted 
to allow for the consequences of Bose-Einstein statistics 
when evaluating the rates of reactions involving multiply- 
deuterated species. 
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H + H -» H 2 (p) 


0.25 






H + H -> H 2 (o) 


0.75 






H + D HD 


1.00 






H + crp — > H+ + e" 


0.46 






Ho + crp — > Hc + + e~ 


0.50 






H 2 (p) + crp -> H+ + H + e- 


0.04 






H 2 (o) + crp -» H+ + H + e- 


0.04 






H 2 (p) + crp -» H + H 


1.50 






H 2 (o) + crp H + H 


1.50 






H 2 (p) + crp -» H 2 + (p) + e- 


0.96 






H 2 (o) + crp -> H 2 +(o) + e" 


0.96 






H+ + H 2 (o) -> H+ + H 2 (p) 


2.20(-10) 


0.00 


0.0 


H 3 +( P ) + H 2 (o) - H 3 +(o) + H 2 (p) 


4.40(-10) 


0.00 


0.0 


H 3 +(o) + H 2 (o) - H 3 +(p) + H 2 (p) 


1.10(-10) 


0.00 


0.0 


H+ + H 2 (p) -f H+ + H 2 (o) 


1.98(-09) 


0.00 


170.5 


H 3 +(p) + H 2 (p) - H 3 +(o) + H 2 (o) 


1.98(-09) 


0.00 


203.4 


H 3 +(o) + H 2 (p) - H 3 +(p) + H 2 (o) 


1.98(-09) 


0.00 


137.6 


H 2 + (o) + H 2 (o) - H 3 +(p) + H 


1.05(-09) 


0.00 


0.0 


H 2 + (o) + H 2 (o) -> H 3 +(o) + H 


1.05(-09) 


0.00 


0.0 


H 2 +(p) + H 2 (o) - H 3 +(p) + H 


1.05(-09) 


0.00 


0.0 


H 2 +(p) + H 2 (o) - H 3 +(o) +H 


1.05(-09) 


0.00 


0.0 


H 2 + (o) +H 2 (p) ^H 3 +(p) +H 


1.05(-09) 


0.00 


0.0 


H 2 +(o) + H 2 (p) -> H 3 +(o) + H 


1.05(-09) 


0.00 


0.0 


H 2 + (p) +H 2 (p) ^H 3 +(p) +H 


2.10(-09) 


0.00 


0.0 


Hc+ + H 2 (p) — > H+ + H + He 


1.10(-13) 


-0.24 


0.0 


He+ + H 2 (o) -> H+ + H + He 


1.10(-13) 


-0.24 


0.0 


H+ + e~ -> H + photon 


3.61(-12) 


-0.75 


0.0 


H 2 + (p) + e- -» H + H 


1.60(-08) 


-0.43 


0.0 


H 2 + (o) + c- -» H + H 


1.60(-08) 


-0.43 


0.0 


He+ + e~ — > He + photon 


4.50(-12) 


-0.67 


0.0 


H 3 +(p) +e ^H + H + H 


5.10(-08) 


-0.52 


0.0 


H 3 +(o) +e ^H + H + H 


5.10(-08) 


-0.52 


0.0 


H 3 +(p) + e- -> H 2 (p) + H 


1.13(-08) 


-0.52 


0.0 


H 3 +(p) + e- -> H 2 (o) + H 


0.57(-08) 


-0.52 


0.0 


H 3 +(o) + e- H 2 (o) + H 


1.70(-08) 


-0.52 


0.0 


D + crp -> D+ + e" 


0.46 






HD + crp -> H+ + D + e" 


0.02 






HD + crp -> D+ + H + e" 


0.02 






HD + crp -> H + D 


1.50 






HD + crp -> HD+ + e" 


0.96 






D+ + H -» H+ + D 


1.00(-09) 


0.00 


0.0 


H+ + D -> D+ + H 


1.00(-09) 


0.00 


41.0 


D+ + H 2 (p) -> H+ + HD 


2.10(-09) 


0.00 


0.0 


D+ + H 2 (o) H+ + HD 


2.10(-09) 


0.00 


0.0 


H+ + HD -» D+ + H 2 (p) 


1.00(-09) 


0.00 


464.0 


H+ + HD -» D+ + H 2 (o) 


1.00(-09) 


0.00 


634.5 


HD+ + H 2 (p) -» H 2 D+(p) + H 


5.25(-10) 


0.00 


0.0 


HD+ + H 2 (p) -> H 2 D+(o) + H 


5.25(-10) 


0.00 


0.0 


HD+ + H 2 (o) -» H 2 D+(p) + H 


5.25(-10) 


0.00 


0.0 


HD+ + H 2 (o) H 2 D+(o) + H 


5.25(-10) 


0.00 


0.0 


HD+ + H 2 (p) -» H 3 +(p) +D 


1.05(-09) 


0.00 


0.0 


HD+ + H 2 (o) H 3 +(o) + D 


5.25(-09) 


0.00 


0.0 


HD+ + H 2 (o) H 3 +(p) + D 


5.25(-09) 


0.00 


0.0 


H 2 + (p) + HD -» H 2 D+(p) + H 


5.25(-10) 


0.00 


0.0 


H 2 + (p) + HD -» H 2 D+(o) + H 


5.25(-10) 


0.00 


0.0 


H 2 + (o) + HD -> H 2 D+(p) + H 


5.25(-10) 


0.00 


0.0 
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p 


H 2 + (o) + HD -> H 2 D+(o) + H 


5.25(-10) 


0.00 


0.0 


H 2 + (p)+HD^H 3 +(p)+D 


1.05(-09) 


0.00 


0.0 


H 2 + (o) + HD -> H 3 +(p) + D 


5.25(-10) 


0.00 


0.0 


H 2 + (o) + HD ^ H 3 +(o) + D 


5.25(-10) 


0.00 


0.0 


H 3 +(p) + D -» H 2 D+(p) + H 


0.66(-09) 


0.00 


0.0 


H 3 + (p) + D -» H 2 D+(o) + H 


0.33(-09) 


0.00 


0.0 


H 3 + (o) + D -» H 2 D+(o) + H 


1.00(-09) 


0.00 


0.0 


H 2 D+(p) + H -> H 3 +(p) + D 


1.00(-09) 


0.00 


632.0 


H 2 D+(o) + H -» H 3 +(o) + D 


0.50(-09) 


0.00 


578.4 


H 2 D+(o) + H -» H 3 +(p) + D 


0.50(-09) 


0.00 


545.5 


H 3 +(p) + HD - H 2 D+(p) + H 2 (p) 


1.17(-10) 


0.00 


0.0 


H 3 +(p) + HD -> H 2 D+(p) + H 2 (o) 


1.17(-10) 


0.00 


0.0 


H 3 +(p) + HD - H 2 D+(o) + H 2 (p) 


5.83(-ll) 


0.00 


0.0 


H 3 +(p) + HD - H 2 D+(o) + H 2 (o) 


5.83(-ll) 


0.00 


25.0 


H 3 +(o) + HD - H 2 D+(p) + H 2 (o) 


1.75(-10) 


0.00 


0.0 


H 3 +(o) + HD - H 2 D+(o) + H 2 (o) 


1.75(-10) 


0.00 


0.0 


H 2 D+(p) + H 2 (o) -» H 2 D+(o) + H 2 (p) 


1.98(-09) 


0.00 


0.0 


H 2 D+(o) + H 2 (p) -> H 2 D+(p) + H 2 (o) 


1.98(-09) 


0.00 


84.0 


H 2 D+(o) + H 2 (o) - H 2 D+(p) + H 2 (p) 


2.44(-ll) 


0.00 


0.0 


H 2 D+(p) + H 2 (p) -» H 2 D+(o) + H 2 (o) 


1.98(-09) 


0.00 


257.0 


H 2 D+(p) + H 2 (p) - H 3 +(p) + HD 


1.40(-10) 


0.00 


232.0 


H 2 D+(p) + H 2 (o) - H 3 +(p) + HD 


0.70(-10) 


0.00 


61.5 


H 2 D+(p) + H 2 (o) -> H 3 +(o) + HD 


0.70(-10) 


0.00 


94.4 


H 2 D+(o) + H 2 (p) -» H 3 +(p) + HD 


1.40(-10) 


0.00 


145.5 


H 2 D+(o) + H 2 (o) - H 3 +(p) + HD 


0.70(-10) 


0.00 


0.0 


H 2 D+(o) + H 2 (o) ^ H 3 +(o) + HD 


0.70(-10) 


0.00 


8.0 


Ho+ + HD -> H+ + D + He 


5.50(-14) 


-0.24 


0.0 


He+ + HD -> D+ + H + He 


5.50(-14) 


-0.24 


0.0 


D+ + e~ -> D + photon 


3.61(-12) 


-0.75 


0.0 


HD+ + e" -» H + D 


3.40(-09) 


-0.40 


0.0 


H 2 D+(p) +e ^H + H + D 


4.96(-08) 


-0.52 


0.0 


H 2 D+(o) +e ^H + H + D 


4.96(-08) 


-0.52 


0.0 


H 2 D+(p) + e - -» HD + H 


1.36(-08) 


-0.52 


0.0 


H 2 D+(o) + e - -> HD + H 


1.36(-08) 


-0.52 


0.0 


H 2 D+(p) + e" -» H 2 (p) + D 


4.76(-09) 


-0.52 


0.0 


H 2 D+(o) + e" -» H 2 (o) + D 


4.76(-09) 


-0.52 


0.0 


D+ + H -> HD+ + photon 


3.90(-19) 


1.80 


0.0 


H+ + D -> HD+ + photon 


3.90(-19) 


1.80 


0.0 


H 2 + (p) + D -> H 2 D+(p) + photon 


7.00(-18) 


1.80 


0.0 


HD+ + H -» H 2 D+(p) + photon 


1.20(-17) 


1.80 


0.0 


HD+ + H -» H+ + HD 


6.40(-10) 


0.00 


0.0 


H 2 + (p) + D - D+ + H 2 (p) 


6.40(-10) 


0.00 


0.0 


HD+ + H -y H 2 +(p) + D 


1.00(-09) 


0.00 


154.0 


H 2 + (p) + D -> HD+ + H 


1.00(-09) 


0.00 


0.0 


D 2 + crp — > D+ + D + e~ 


0.04 






D 2 + crp -> D + D 


1.50 






D 2 + crp -> D 2 + + e" 


0.96 






H+ + D 2 -> D+ + HD 


2.10(-09) 


0.00 


491.0 


D+ + HD -> H+ + D 2 


1.00(-09) 


0.00 


0.0 


D 2 + + H 2 (p) -> H 2 D+(p) + D 


1.05(-09) 


0.00 


0.0 


D 2 + + H 2 (o) -» H 2 D+(o) + D 


1.05(-09) 


0.00 


0.0 


D 2 + + H 2 (p) -> HD 2 + + H 


1.05(-09) 


0.00 


0.0 


D 2 + + H 2 (o) -» HD 2 + + H 


1.05(-09) 


0.00 


0.0 


H 2 + (p) + D 2 -> H 2 D+(p) + D 


1.05(-09) 


0.00 


0.0 


H 2 + (o) + D 2 -» H 2 D+(o) + D 


1.05(-09) 


0.00 


0.0 


H 2 + (p) + D 2 -> HD 2 + + H 


1.05(-09) 


0.00 


0.0 
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a 


H 2 + (o) + D 2 -» HD 2 + + H 


1.05( 


-09) 


0.00 


0.0 


HD+ + HD -> H 2 D+(p) + D 


5.25( 


-10) 


0.00 


0.0 


HD+ + HD -> H 2 D+(o) + D 


5.25( 


-10) 


0.00 


0.0 


HD+ + HD -> HD 2 + + H 


1.05( 


-09) 


0.00 


0.0 


H 2 D+(p) + D -> HD 2 + + H 


1.00( 


-09) 


0.00 


0.0 


H 2 D+(o) + D -f HD 2 + + H 


1.00( 


-09) 


0.00 


0.0 


HD 2 + + H -> H 2 D+(p) + D 


1.00( 


-09) 


0.00 


600.0 


HD 2 + + H -> H 2 D+(o) + D 


1.00( 


-09) 


0.00 


513.5 


H 3 +(o) + D 2 -f H 2 D+(o) + HD 


3.50( 


-10) 


0.00 


0.0 


H 3 +(o) + D 2 ^ HD 2 + + H 2 (o) 


3.50( 


-10) 


0.00 


0.0 


H 3 +( P ) + D 2 - HD 2 + + H 2 (p) 


2.33( 


-10) 


0.00 


0.0 


H 3 +(p) + D 2 -> HD 2 + + H 2 (o) 


1.17( 


-10) 


0.00 


0.0 


H 3 +(p) + D 2 -» H 2 D+(p) + HD 


2.33( 


-10) 


0.00 


0.0 


H 3 +(p) + D 2 -> H 2 D+(o) + HD 


1.17( 


-10) 


0.00 


0.0 


H 2 D+(p) + HD ^ HD 2 + + H 2 (p) 


2.60( 


-10) 


0.00 


0.0 


H 2 D+(o) + HD -> HD 2 + + H 2 (o) 


2.60( 


-10) 


0.00 


0.0 


H 2 D+(o) + HD -> H 3 + (o) + D 2 


1.30( 


-10) 


0.00 


108.0 


H 2 D+(o) + HD -» H 3 +(p) + D 2 


1.30( 


-10) 


0.00 


75.0 


H 2 D+(p) + HD -» H 3 +(p) + D 2 


2.60( 


-10) 


0.00 


161.5 


HD 2 + + H 2 (p) -» H 2 D+(p) + HD 


1.00( 


-10) 


0.00 


190.0 


HD 2 + + H 2 (p) -» H 2 D+(o) + HD 


1.00( 


-10) 


0.00 


277.0 


HD 2 + + H 2 (o) H 2 D+(o) + HD 


1.00( 


-10) 


0.00 


106.0 


HD 2 + + H 2 (o) ^ H 2 D+(p) + HD 


1.00( 


-10) 


0.00 


20.1 


HD 2 + + H 2 (p) -» H 3 +(p) + D 2 


2.00( 


-10) 


0.00 


320.0 


HD 2 + + H 2 (o) H 3 +(p) + D 2 


1.00( 


-10) 


0.00 


149.5 


HD 2 + + H 2 (o) H 3 +(o) + D 2 


1.00( 


-10) 


0.00 


182.4 


Hc+ + D 2 D+ + D + Ho 


1.10( 


-13) 


-0.24 


0.0 


Hc+ + D 2 D 2 + + He 


2.50( 


-14) 


0.00 


0.0 


D 2 + + e~ -> D + D 


3.40( 


-09) 


-0.40 


0.0 


HD 2 + + e -->D + D + H 


4.96( 


-08) 


-0.52 


0.0 


HD 2 + + e " -> HD + D 


1.36( 


-08) 


-0.52 


0.0 


HD 2 + + e - -> D 2 + H 


4.76( 


-09) 


-0.52 


0.0 


D 2 + + H -> H+ + D 2 


6.40( 


-10) 


0.00 


0.0 


HD+ + D -> D+ + HD 


6.40( 


-10) 


0.00 


0.0 


HD+ + D -> D 2 + + H 


1.00( 


-09) 


0.00 


0.0 


D 2 + + H -> HD+ + D 


1.00( 


-09) 


0.00 


472.0 


D 2 + + D D+ + D 2 


6.40( 


-10) 


0.00 


0.0 


D 2 + + HD HD 2 + + D 


1.05( 


-09) 


0.00 


0.0 


D 2 + + HD D 3 + + H 


1.05( 


-09) 


0.00 


0.0 


HD+ + D 2 HD 2 + + D 


1.05( 


-09) 


0.00 


0.0 


HD+ + D 2 D 3 + + H 


1.05( 


-09) 


0.00 


0.0 


HD 2 + + D — ► D 3 + + H 


1.00( 


-09) 


0.00 


0.0 


D 3 + + H -> HD 2 + + D 


1.00( 


-09) 


0.00 


655.0 


HD 2 + + HD ^ D 3 + + H 2 (p) 


1.00( 


-10) 


0.00 


0.0 


HD 2 + + HD -> D 3 + + H 2 (o) 


1.00( 


-10) 


0.00 


0.0 


HD 2 + + HD H 2 D+(p) + D 2 


1.00( 


-10) 


0.00 


110.0 


HD 2 + + HD -► H 2 D+(o) + D 2 


1.00( 


-10) 


0.00 


196.5 


H 2 D+(p) + D 2 HD 2 + + HD 


8.50( 


-10) 


0.00 


0.0 


H 2 D+(o) + D 2 HD 2 + + HD 


8.50( 


-10) 


0.00 


0.0 


H 2 D+(p) + D 2 - D 3 + + H 2 (p) 


8.50( 


-10) 


0.00 


0.0 


H 2 D+(o) + D 2 -► D 3 + + H 2 (o) 


8.50( 


-10) 


0.00 


0.0 


D 3 + + H 2 (p) -» H 2 D+(p) + D 2 


1.50( 


-09) 


0.00 


352.0 


D 3 + + H 2 (o) H 2 D+(o) + D 2 


1.50( 


-09) 


0.00 


268.0 


D 3 + + H 2 (p) -» HD 2 + + HD 


1.50( 


-09) 


0.00 


241.0 


D 3 + + H 2 (o) HD 2 + + HD 


1.50( 


-09) 


0.00 


70.5 


D 3 + + e~^D + D + D 


2.03( 


-08) 


-0.50 


0.0 
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7 


a 


a 


D 3 + + e- -> D 2 + D 


6.75( 


-09) 


-0.50 


0.0 


D 2 + + D 2 -> D 3 + + D 


2.10( 


-09) 


0.00 


0.0 


D 3 + + HD -> HD 2 + + D 2 


1.50( 


-09) 


0.00 


161.5 


g° + sccpho — > g+ + e~ 


0.63 


08) 






g~ + sccpho — > g° + e~ 


0.41 


09) 






g° + e~ — > g~ + photon 


6.90( 


-05) 


0.50 


0.0 


g- + H+ - g° + H 


1.60( 


-06) 


0.50 


0.0 


g + H 3 +(p) -> g° + H 2 (p) + H 


3.07( 


-07) 


0.50 


0.0 


g + H 3 +(p) - g° + H 2 (o) + H 


1.54( 


-07) 


0.50 


0.0 


g- + H 3 +(o) -> g° + H 2 (o) + H 


4.61( 


-07) 


0.50 


0.0 


g- + H 3 +(p) -> g° + 3H 


4.61( 


-07) 


0.50 


0.0 


g- + H 3 +(o) - g° + 3H 


4.61( 


-07) 


0.50 


0.0 


g" + Hc+ -> g° + He 


8.00( 


-07) 


0.50 


0.0 


g° + H+ -> g+ + H 


1.60( 


-06) 


0.50 


0.0 


g° + H 3 +(p) -> g+ + H 2 (p) + H 


3.07( 


-07) 


0.50 


0.0 


g° + H 3 +(p) -> g+ + H 2 (o) H 


1.54( 


-07) 


0.50 


0.0 


g° +H 3 +(o) ^g+ + H 2 (o) H 


4.61( 


-07) 


0.50 


0.0 


g° + H 3 +(p) -> g+ + 3H 


4.61( 


-07) 


0.50 


0.0 


g° + H 3 +(o) -> g+ + 3H 


4.61( 


-07) 


0.50 


0.0 


g° + Hc+ -> g+ + He 


8.00( 


-07) 


0.50 


0.0 


g + + c~ — > g° + photon 


6.90( 


-05) 


0.50 


0.0 


g- + D+ -> g° + D 


1.13( 


-06) 


0.50 


0.0 


g- + H 2 D+(p) ^ g° + H 2 (p) + D 


1.33( 


-07) 


0.50 


0.0 


g + H 2 D+(p) -» g° + HD + H 


2.66( 


-07) 


0.50 


0.0 


g- + H 2 D+(o) ^ g° + H 2 (o) + D 


1.33( 


-07) 


0.50 


0.0 


g + H 2 D+(o) ^ g° + HD + H 


2.66( 


-07) 


0.50 


0.0 


g- + H 2 D+(p) -» g° + D + 2H 


3.99( 


-07) 


0.50 


0.0 


g- + H 2 D+(o) -> g° + D + 2H 


3.99( 


-07) 


0.50 


0.0 


g° + D+ -» g+ + D 


1.13( 


-06) 


0.50 


0.0 


g° + H 2 D+(p)^g+ + H 2 (p) + D 


1.33( 


-07) 


0.50 


0.0 


g° + H 2 D+(p) -» g+ + HD + H 


2.66( 


-07) 


0.50 


0.0 


g° + H 2 D+(o) ^ g+ + H 2 (o) + D 


1.33( 


-07) 


0.50 


0.0 


g° + H 2 D+(o) -» g+ + HD + H 


2.66( 


-07) 


0.50 


0.0 


g° + H 2 D+(p) -» g+ + D + 2H 


3.99( 


-07) 


0.50 


0.0 


g° + H 2 D+(o) -» g+ + D + 2H 


3.99( 


-07) 


0.50 


0.0 


g + HD 2 + ^ g° + HD + D 


2.38( 


-07) 


0.50 


0.0 


g- + HD 2 + ^ g° + D 2 + H 


1.19( 


-07) 


0.50 


0.0 


g + HD 2 + ^ g° + 2D + H 


3.57( 


-07) 


0.50 


0.0 


g° + HD 2 + -» g+ + HD + D 


2.38( 


-07) 


0.50 


0.0 


g° + HD 2 + -» g+ + D 2 + H 


1.19( 


-07) 


0.50 


0.0 


g° + HD 2 + -> g+ + 2D + H 


3.57( 


-07) 


0.50 


0.0 


g- + D 3 + -» g° + D 2 + D 


3.26( 


-07) 


0.50 


0.0 


g- + D 3 + -» g° + 3D 


3.26( 


-07) 


0.50 


0.0 


g° + D 3 + -> g+ + D 2 + D 


3.26( 


-07) 


0.50 


0.0 


g° + D 3 + g+ + 3D 


3.26( 


-07) 


0.50 


0.0 
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Table A.l: Reactions and rate coefficients adopted in our chem- 
ical model. The parameters a, (3, and 7 define the rate coef- 
ficients k cm 3 s _1 at temperature T through the relation k = 
3(a g , T)7(T/300) a exp(-/3/T); J allows for Coulomb focusing in re- 
actions of positive ions and negatively charged grains (Draine and 
Sutin 1987, equation 3.4); a g = 0.1 /im is adopted in the Table. The 
rates s _1 of reactions induced directly by cosmic rays (crp) or indi- 
rectly by the secondary electrons, which generate H2 fluorescence 
photons (secpho), are given by 7C, where £ is the rate of cosmic 
ray ionization of H2. In the first three reactions (formation of H2 
or HD on grain surfaces), 7 denotes the fraction of such reactions 
which form the specified product; the rate is calculated internally 
in the program from the grain parameters. Key reactions arc in 
bold face. Numbers in parentheses are powers of 10. 
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Appendix B: Ice Mantle Thickness 

In this Appendix, we derive estimates of the ice-mantle 
thickness in completely depleted cores. We follow the ap- 
proach of FPdesF. 

Following FPdesF, we assume the refractory cores to 
be composed of silicates (in practice, olivine, MgFeSiC^) 
or graphite. The mass of the cores, per unit volume of 
gas, is n{X)M, where the summation extends over the 
elements X, and M denotes the elemental mass. We define 
a c through the relation n(X)M = a c bhihh, and a c = 
7.84 10 -3 , which corresponds to a fraction of 0.0056 by 
mass relative to the gas-phase H and He. 

Assuming that all elements heavier than He are in the 
grain cores or frozen in the mantles, the total dust-to- 
gas mass density ratio is 0.0125, and the mass per unit 
volume of the gas of the ice in the mantle, is cvjnHTtH, 
where a, = 9.63 10~ 3 . 

The volume V~i of a grain ice mantle is determined by 
the relation 

n c ViPi = a^HmH (B.l) 

where n c is the number of grain cores (i.e. the number of 
grains) per unit volume of the gas and pi is the density of 
the mantle material; n c is defined by 

n c V c p c = a c n H m H (B.2) 

where p c is the density and V c is the volume of refractory 
core material. We adopt p c = 3 gm cm -3 and pi — 1 gm 
cm~ 3 . Hence, the ratio of mantle to core volume is 

Vt /V c = ^ (B.3) 

OtcPi 

Substituting the above numbers, we find that the ra- 
tio VijV c = 3.72, or a fraction 0.79 of the grain volume is 
occupied by the mantle. The ratio of mantle thickness to 
core radius is 0.67, and the ratio of mantle thickness to 
grain (core + mantle) radius is 0.40, in the limit of com- 
plete heavy element depletion. Thus, for a grain of radius 
a g = 0.1 /um, the mantle thickness is 400 A. 

This analysis can be generalized to the case of a MRN 
power-law size distribution, using the formalism of Le 
Bourlot et al. (1995). If the lower limit of the size dis- 
tribution is taken to be a m = 0.01 ^m (and the upper 
limit is a,M = 0.3 /nm), then the thickness 6 of the mantle 
is obtained by solving the equation 

U{Vi/V c ) = (/ 1( 5 3 + 3/ 2 <5 2 + 3/ 3 <f) (B.4) 

where /„ is defined by 

/„ = (a n M r a^ r )/(n r) (B.5) 

and r = 3.5 for the MRN grain size distribution. Solving 
for 5, we find that the mantle thickness is 260 A, indepen- 
dent of a for the specified values of a m (and Om)- Thus, 
very small grains are not present under conditions of com- 
plete heavy element depletion, owing to the formation of 
an ice mantle around the core. 
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